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^ • ABSTRACT 

i^'. Galaxy clustering and galaxy-galaxy lensing probe the connection between galaxies 

^ . and their dark matter haloes in complementary ways. Since the clustering of dark 

H ' matter haloes depends on cosmology, the halo occupation statistics inferred from the 

C/3 , observed clustering properties of galaxies are degenerate with the adopted cosmology. 

^ ^^ ■ Consequently, different cosmologies imply different mass-to-light ratios for dark matter 

haloes. Galaxy-galaxy lensing, which yields direct constraints on the actual mass-to- 
light ratios, can therefore be used to break this degeneracy, and thus to constrain 
^ ■ cosmological parameters. In this paper we establish the link between galaxy lumi- 

0\J I nosity and dark matter halo mass using the conditional luminosity function (CLF), 

CO i 4>(L|M )dL, which gives the number of galaxies with luminosities in the range L±dL/2 

0^ ' that reside in a halo of mass M. We constrain the GLF parameters using the galaxy 

\l^ I luminosity function and the luminosity dependence of the correlation lengths of galax- 

[*-. ■ ies. The resulting CLF models are used to predict the galaxy-galaxy lensing signal. 

f^ ' For a cosmology that agrees with constraints from the cosmic microwave background, 

OO , i.e. (rim, CTs) = (0.238, 0.734), the model accurately fits the galaxy-galaxy lensing data 

^D ■ obtained from the SDSS. For a comparison cosmology with (f^m, as) = (0-3, 0.9), how- 

ever, we can accurately fit the luminosity function and clustering properties of the 
galaxy population, but the model predicts mass-to-light ratios that are too high, re- 
/\^ ' suiting in a strong overprediction of the galaxy-galaxy lensing signal. We conclude that 

JH I the combination of galaxy clustering and galaxy-galaxy lensing is a powerful probe of 

the galaxy-dark matter connection, with the potential to yield tight constraints on 
cosmological parameters. Since this method mainly probes the mass distribution on 
relatively small (non-linear) scales, it is complementary to constraints obtained from 
the galaxy power-spectrum, which mainly probes the large-scale (linear) matter dis- 
tribution. 

Key words: galaxies: halos — large-scale structure of Universe — dark matter — 
cosmological parameters — gravitational lensing — methods: statistical 



1 INTRODUCTION Since galaxies are believed to form and reside in dark mat- 

ter haloes, the clustering strength of a given population of 
With the advent oi large galaxy redshnt surveys, it has be- , . , i . . i . r i i . . i i 

galaxies can be compared to that ot dark matter haloes as 
come possible to obtain accurate measurements ot the cms- i- j_ i i ■ i ■ i j_- ^i j_ i i n 

predicted by numerical simulations or the extended Press- 
tering of galaxies as a function of their properties, such as In- o i i ^ r i- o i ■ i ^^^ r 

" ° , , / ^ ^T achechter lormalism. buch a comparison reveals a wealth ot 

mmosity, morphology and color (e.g. Guzzo et al. 2000; Nor- . „ ,. , , ,, hit 77 

^ , , „, , , mtormation about the so-called qalaxy-dark matter connec- 

berg et al. 2001, 2002; Zehavi et al. 2005; Wang et al. 2007 . ^. . t- a^ p -o- mno 1 p q -^u onnn 

^ ' 10 / jjQ^ (-g g jing, Mo & corner 1998; Peacock & Smith 2000; 

Berlind & Weinberg 2002; Yang, Mo & van den Bosch 2003; 

van den Bosch, Yang & Mo 2003; van den Bosch et al. 2007). 
* International Max-Planck Research School Fellow Unfortunately, this method of constraining the link be- 
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tering has one important shortcoming: the halo occupation 
statistics inferred from the observed clustering properties 
depend on the cosmological parameters adopted. More pre- 
cisely, models based on different cosmologies can fit the clus- 
tering data equally well by simply relying on different halo 
occupation statistics or, equivalently, different mass-to-light 
ratios. In order to break this degeneracy between cosmol- 
ogy and halo occupation statistics independent constraints 
on the mass-to-light ratios are required (e.g. van den Bosch, 
Mo & Yang 2003; Tinker et al. 2005). One method that 
can provide these constraints is galaxy-galaxy lensing (here- 
after g-g lensing), which probes the mass distributions (and 
hence the halo masses) around galaxies. This implies that 
the combination of clustering and lensing in principle holds 
the potential to put constraints on cosmological parameters 
(Seljak et al. 2005; Yoo et al. 2006). 

The first attempt to detect g-g lensing was made by 
Tyson et al. (1984), but because of the relatively poor qual- 
ity of their data they were unable to detect a statistically sig- 
nificant signal. With the advent of wider and deeper surveys 
becoming available, however, g-g lensing has now been de- 
tected with very high significance, and as function of various 
properties of the lensing galaxies (e.g. Griffiths et al. 1996 
Hudson et al. 1998; McKay et al. 2001; Guzik & Seljak 2002 
Hoekstra et al. 2003, 2004; Sheldon et al. 2004, 2007a,b: 
Mandelbaum et al. 2006; Heymans et al. 2006; Johnston 
et al. 2007; Parker et al. 2007; Mandelbaum, Seljak & Hi- 
rata 2008). Unfortunately, a proper interpretation of these 
data in terms of the link between galaxies and dark mat- 
ter haloes has been hampered by the fact that the lensing 
signal can typically only be detected when stacking the sig- 
nal of many lenses. Since not all lenses reside in haloes of 
the same mass, the resulting signal is a non-trivial aver- 
age of the lensing signal due to haloes of different masses. 
Most studies to date have assumed that the relation between 
the luminosity of a lens galaxy and the mass of its halo is 
given by a simple power-law relation with zero scatter (see 
Limousin et al. 2007 for a detailed overview). However, it 
has become clear, recently, that the scatter in this relation 
between light and mass can be very substantial (More et 
al. 2008b, and references therein). As shown by Tasitsiomi 
et al. (2004), this scatter has a very significant impact on 
the actual lensing signal, and thus has to be accounted for 
in the analysis. In addition, central galaxies (those residing 
at the center of a dark matter halo) and satellite galaxies 
(those orbiting around a central galaxy) contribute very dif- 
ferent lensing signals, even when they reside in haloes of 
the same mass (e.g. Natarajan, Kneib & Small 2002; Yang 
et al. 2006; Limousin et al. 2007). This has to be properly 
accounted for (see e.g. Guzik & Seljak 2002), and requires 
knowledge of both the satellite fractions and of the spatial 
number density distribution of satellite galaxies within their 
dark matter haloes. 

Over the years, numerous techniques have been de- 
veloped to interpret galaxy-galaxy lensing measurements 
(Schneider & Rix 1997; Natarajan & Kneib 1997; Brain- 
erd & Wright 2002; Guzik & Seljak 2001). Several authors 
have also used numerical simulations to investigate the link 
between g-g lensing and the galaxy-dark matter connnection 
(e.g., Tasitsiomi et al. 2004; Limousin et al. 2005; Natarajan, 
De Lucia & Springel 2007; Hayashi & White 2007). It has 
become clear from these studies that g-g lensing in principle 



contains a wealth of information regarding the mass distri- 
butions around galaxies; in addition to simply probing halo 
masses, g-g lensing also holds the potential to measure the 
shapes, concentrations and radii of dark matter haloes, and 
the first observational results along these lines have already 
been obtained (Natarajan et al. 2002; Hoekstra et al. 2004; 
Mandelbaum et al. 2006; Limousin et al. 2007; Mandelbaum 
et al. 2008). 

In this paper we use an analytical model, similar to that 
developed by Seljak (2000) and Guzik & Seljak (2001), to 
predict the g-g lensing signal as a function of the luminosity 
of the lenses starting from a model for the halo occupation 
statistics that is constrained to fit the abundances and clus- 
tering properties of the lens galaxies. A comparison of these 
predictions with the data thus allows us to test the mass- 
to-light ratios inferred from the halo occupation model, and 
ultimately to constrain cosmological parameters. The model 
assumes that haloes have NFW (Navarro, Frenk & White 
1997) density distributions, and that satellite galaxies fol- 
low a radial number density distribution that is unbiased 
with respect to the dark matter. The occupation statistics 
are described via the conditional luminosity function (CLF; 
see Yang et al. 2003) , which specifies the average number of 
galaxies of given luminosity that reside in a halo of given 
mass. This CLF is ideally suited to model g-g lensing, as it 
allows one to properly account for the scatter in the relation 
between luminosity and halo mass, and to split the galaxy 
population in centrals and satellites. We demonstrate how 
these different galaxy populations contribute to the lensing 
signal in different luminosity bins, and show that uncertain- 
ties related to the expected concentrations of dark matter 
haloes and the radial number density distributions of satel- 
lite galaxies do not have a significant impact on our results. 

Assuming a fiat ACDM cosmology with parameters sup- 
ported by the third year data release of the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP, see Spergel et al. 2007), 
we obtain a CLF that accurately fits the abundances and 
clustering properties of SDSS galaxies. Using our analytical 
model, we show that this same CLF also accurately matches 
the g-g lensing data obtained from the SDSS by Seljak et 
al. (2005) and Mandelbaum et al. (2006) without any addi- 
tional tuning of the model parameters. However, if we repeat 
the same exercise for a cosmology with a matter density and 
power-spectrum normalization that are slightly (~ 20 per- 
cent) higher, the model that fits the clustering data can 
no longer simultaneously fit the g-g lensing data. This con- 
firms that a joint analysis of clustering and g-g lensing is an 
extremely promising method to constrain cosmological pa- 
rameters. In a companion paper (Li et al. 2008), we use the 
SDSS galaxy group catalogue of Yang et al. (2007) to pre- 
dict the g-g lensing signal, which we compare to data from 
the SDSS. Although, Li et al. obtain their halo occupation 
statistics from a galaxy group catalogue, rather than from 
the galaxy clustering properties, they obtain very similar 
results. 

The present paper is organized as follows. We review 
the necessary formalism of g-g lensing in §[3 with a detailed 
description of the model used to interpret the g-g lensing 
signal. The CLF, used to describe the connection between 
galaxies and dark matter haloes, is introduced in § O The 
properties of the predicted g-g lensing signal are illustrated 
in § [4] together with a comparison between theoretical pre- 



Clustering and Galaxy- Galaxy Lensing 3 



dictions and SDSS data. A detailed analysis of the assump- 
tions entering the model is presented in § [S] Conclusions are 
presented in § [6l 



2 THE HALO-MODEL DESCRIPTION OF 
GALAXY-GALAXY LENSING 

Galaxy-galaxy lensing measures the tangential shear distor- 
tions, 7t, in the shapes of background galaxies (hereafter 
sources) induced by the mass distribution around foreground 
galaxies (hereafter lenses). Since the tangential shear dis- 
tortions due to a typical lens galaxy (and its associated 
dark matter halo) are extremely small, and since back- 
ground sources have non-zero intrinsic ellipticities, measur- 
ing 7t with sufficient signal-to-noise requires large numbers 
of background galaxies. Except for extremely deep surveys 
behind clusters of galaxies, which have a large surface area, 
the number density of background sources is insufficient for a 
reliable measurement of 7t around individual lenses. In prac- 
tice, this problem can be circumvented by stacking many 
lenses according to some observable property. For exam- 
ple, Mandelbaum et al. (2006) measured 7t as a function of 
the transverse comoving distance R by stacking thousands 
of lenses in a given luminosity bin [Li,L2]. The resulting 
shear ■yt{R\Li,L2) holds information regarding the charac- 
teristic mass of the haloes that host galaxies with luminos- 
ity Li ^1/^1/2, and hence can be used to constrain the 
galaxy-dark matter connection. 

The tangential shear as a function of the projected ra- 
dius R around the lenses is related to the excess surface 
density (ESD) profile, AE(_R), according to 



AE(i?) = E(< R) - J:{R) = 7t(i?)Sc 



(1) 



where E(_R) is the projected surface density and E(< R) is 
its average inside R, 



-)/\ 7->/ T 7-)/ 



E(< R) = — I E{R')R'dR 



(2) 



(Miralda-Escude 1991; Sheldon et al. 2004). The so-called 
critical surface density, Ecrit, is a geometrical parameter 
given by 

„2 



Ecrit 



WS 



4:TVG OJLOJLsi^ + Zl) ' 



(3) 



with ws, i^L and ojls the comoving distances to the source, 
the lens and between the two, respectively, and with z^ the 
redshift of the lens. 

The projected surface density is related to the galaxy- 
dark matter cross correlation, ^g,dm(''), according to 



E(i?) 



[l + ^g,dm(r-)] Au, 



(4) 



where p is the average density of matter in the Universe and 
the integral is along the line of sight with uj the comoving 
distance from the observer. The three-dimensional comoving 
distance r is related to a; through r^ = luI^+uj'^ — 2uj-luj cos d 
(see Fig.[l]for an illustration of the geometry). Since ^g,dm(7') 
goes to zero in the limit r ^ oo, and since in practice 6 is 
small, we can approximate Eq. (|4]) using 



E(i?) = 2p/ [1 + Cg,dm(r)] 



rdr 



Vr2 - R2 



(5) 




Figure 1. Illustration of the geometry between source, lens and 
observer 



which is the expression we adopt throughout. 

The main goal of this paper is to test the halo occu- 
pation statistics inferred from galaxy clustering data with 
g-g lensing data. As is evident from the above equations, 
the lensing signal AE(_R) is completely specified by the 
galaxy dark matter cross correlation, which, as we demon- 
strate below, can be computed from a given halo occupation 
model. For computational convenience, we will be working in 
Fourier space, where the related quantity is the galaxy-dark 
matter cross power spectrum 



-Pg,dm(fc) =47r / Cg,dm('-) 



f 

Jo 



sin(fcr) 



kr 



r dr . 



(6) 



In order to compute this power spectrum, we follow Seljak 
(2000) and Guzik & Seljak (2001), and adopt the halo model, 
according to which all dark matter is partitioned over dark 
matter haloes (see also Mandelbaum et al. 2005a) . As usual 
in the halo model, it is convenient to split Pg,dm(fe) into two 
terms; a 1-halo term, which describes the cross correlation 
between galaxies and the dark matter particles that reside 
in the same halo, and a 2-halo term, where each galaxy is 
cross correlated with the dark matter in all haloes except 
for the one that hosts the galaxy in question. The compu- 
tation of these two terms has to account for two important 
complications. First of all, because of the stacking procedure 
used in order to achieve sufficient signal-to-noise, the ESD 
contains signal from haloes with different masses. A proper 
estimate of Pg,dm(fc), therefore, requires the full probability 
distribution that a galaxy with the stacking property used 
(in this case luminosity) resides in a dark matter halo of 
mass M. Secondly, central galaxies (those residing at the 
center of a dark matter halo) and satellite galaxies (those 
orbiting around a central galaxy) contribute very different 
lensing signals, even when they reside in haloes of the same 
mass (e.g., Yang et al. 2006). This has to be properly ac- 
counted for, and requires knowledge of both the satellite 
fractions and of the spatial number density distribution of 
satellite galaxies within their dark matter haloes. Based on 
these considerations, we split Pg,dm(fc) in four terms: 



\P',X(k) 



p!yjk) 



+ .h[PX(k)+P!X(k)] 



(7) 



where 'c' and 's' stand for 'central' and 'satellite', respec- 
tively. The reason for explicitely writing the central and 
satellite fractions (/c and fa = 1 — /c, respectively) in the 
above equation will become apparent below, in which we 
describe each of these four terms in turn. 
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2.1 The 1-halo term 

The one-halo central term of the power spectrum describes 
the dark matter distribution inside haloes hosting central 
galaxies. For a single, central lensing galaxy, it simply re- 
flects the Fourier transform of the overdensity of the dark 
matter halo in which the lens resides: 



M 



(8) 



where iidm(fc|M) is the normalized Fourier transform of the 
mass density profile, p{r-\M): 

^dm(fc|M) =47r / ^^ ' ' , ' rdr. (9) 



M 



kr 



with ri8o the radius of the halo (see i)2. 31 below). However, 
because the lensing signal is measured by stacking galaxies 
with luminosities in the range [Li, L2], we have that 

^g dm(fc) = = / Vc{M\Li , L2) M ui^{k\M) dM , (10) 
P Jo 

where VdMlLi, L2) is the probability that a central galaxy 
with luminosity Li ^ L ^ L2 resides in a halo of mass 
M. This probability function reflects the halo occupation 
statistics, and, using Bayes' theorem, can be written as 

P.(M|L,,L2)dA/=i^Mll(^ii4M^dM. (11) 

Here {Nc)m{Li,L2) is the average number of central galax- 
ies with luminosities in the range [Li,I/2] that reside in a 
halo of mass M, n{M) is the halo mass function and 



nc(Li,L2)=/ {N,)M{Li,L2)n{M)dM, (12) 

Jo 

is the comoving number density of central galaxies in the 
given luminosity range. 

Combining (|lUp and pip , the first term of the galaxy- 
dark matter power spectrum can be written as 



fcP'X.ik) = -: 



1 



ritotp 

{iVc)M(Li, L2) MMdm(fe|M) n{M) dM (13) 



where ritot ~ nc(Z/i, L2)//c is the total number density of 
all galaxies (centrals plus satellites) with luminosities in the 
range [Li,L2]. Note that, for brevity, we don't explicitely 
write the luminosity dependence of /c and ntot, but it is 
understood that /c = /c(Li, L2) and ntot = ntot{Li,L2). 

The 1-halo satellite term is similar to the 1-halo cen- 
tral term, except for the fact that satellite galaxies do not 
reside at the center of their dark matter halo, but follow 
a number density distribution ns{r\M). Consequently, the 
1-halo lensing signal due to satellite galaxies involves a con- 
volution of ns(r\AI) with the mass density profile p{r\Ad) 
of the host halo in which they reside. Using that in Fourier 
space a convolution corresponds to a simple multiplication, 
we obtain: 

^gdmW = = / r,iM\L^,L2)u,ik\M)Mu^^{k\M)dM, 
P Jo 



with 



Us{k\M) = An 



ns{r\M) 



sin(fcr) ■, , 
^ ^ rdr 



{Ns)m{Li,L2) kr 



(14) 



(15) 



the Fourier transform of ns(r|M) normalized by 
{Ns)m{Li, L2) which is the average number of satel- 
lites with Li ^ L !^ L2 that reside in a halo of mass M. 
We assume that there is no luminosity segregation amongst 
the satellites, so that they all follow the same radial profile, 
independent of their luminosity. We write the probability 
that a satellite galaxy with Li ^ L $; L2 resides in a halo 
of mass M as 



VdM\L,,L2)dM ^ ^-^^^^4^^1l-^^^¥^dM 



ns{Li,L2 



(16) 



with 



ns{Li,L2 



f 

Jo 



{Ns)MiLuL2)n{M)dM, (17) 



the comoving number density of satellite galaxies with lumi- 
nosities in the range [Li,L2]. The 1-halo satellite term can 
thus be written as 



ntotp 



{Ns}m{Li,L2) u,{k\M) M Ui^{k\M) n{M) dM . 

(18) 

where we have used that rTtot = ns(I/i, I/2)//s. Note that 
/, = /,(Li,L2). 

2.2 The 2- halo term 

The 2-halo term of the power spectrum describes the corre- 
lation between galaxies and dark matter particles belonging 
to separate haloes. Within the halo model, this means cross 
correlating each galaxy with all the dark matter haloes other 
than the one in which the galaxy in question resides. Using 
the fact that dark matter haloes are a biased tracer of the 
dark matter mass distribution, the contribution to the 2-halo 
term due to central galaxies can be written as 



g,dm 



(fc) = 



-■dm (k) 



I 



Vc{M\Li,L2)b{M)dM 



P JO 

PI' Mdm(fclM') b{M') n{M') dM' . (19) 



where P^(fc) and b{M) are the non-linear power spectrum 
of the dark matter and the halo bias function, respectively. 
The first integral refiects the contribution of the central 
galaxies, while the second integral describes the dark matter 
density field partitioned over haloes. Using (|lip we obtain 

dNL/ 
dm I 

ntot'p Jo 

M' Udm{k\M') b{M') n(M') dM' . (20) 



/c^gdm«=^^^ / {iVc)M(Li,L2)&(M)n(M)dM 






f^P^Xi^) = ^^^ / ^^' «dm(fc|M') fo(M') niM') dM' 



Similarly, the satellite part of the 2-halo term is given by 

PZjk) f 
ntot'p Jo 

{N^)m{Li,L2) Us{k\M) b{M) n{M) dM . (21) 

where the second integral now accounts for the number den- 
sity distribution of satellite galaxies in haloes of mass M. 

Note that equations (|20p and (|2ip ignore halo exclusion, 
i.e. the fact that, in the halo model, haloes can not overlap. 
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In the Appendix, we present an approximate method to take 
halo exclusion into account. Far from being a detailed treat- 
ment, the suggested procedure accounts only for the most 
relevant effect, i.e. the exclusion of dark matter particles 
residing in the same host halo of central galax;ies (see Ap- 
pendix for further details). Unless stated otherwise, all the 
results shown throughout the paper are obtained applying 
halo exclusion as modelled in the Appendix. 

In addition, a technical, as well as conceptual, issue 
arises in calculating the 2-halo terms introduced in equa- 
tions (|20p and (|2ip . Let us rewrite these two equations in 
the following compact form: 



Table 1. Cosmological Parameters 



g,dm 

p2h,s 

g,dm 



(fc) 
(k) 



= -Prtm(fc)2rjV, 



iMik) 



dm! 



(22) 



where XjVc iN^ik) and lM{k) are 



Jo ^s 

TM{k)^ ^uam{k\M)b{M)n{M)dM . (23) 

Jo P 

The evaluation of these integrals is somewhat tedious nu- 
merically, as it requires knowledge of the halo mass func- 
tion and the halo bias function over the entire mass range 
[0, cxj). Since these have only been tested against nu- 
merical simulations over a limited range of halo masses 
(10^/i"^Mo < M < 10^^ H-'^Mq), it is also unclear how ac- 
curate they are. In practice, though, these problems can be 
circumvented as follows. First of all, because of the exponen- 
tial cut-off in the halo mass function, it is sufficiently accu- 
rate to perform the integrations of Eq. (|23p only up to M — 
10^^h~^ Mq. Secondly, Xjv^ and XN^ik) contain the halo oc- 
cupation statistics, {Nc)m{Li,L2) and (Ns) m (Li , L2) , re- 
spectively, which, for all luminosities of interest in this pa- 
per, are equal to zero for M < 10^h~^ Mq. Therefore, Xn^ 
and XNs{k) can be computed accurately by only integrating 
over the mass range [10^ — 10^^]h~^ Mq. Unfortunately, the 
integrand of TM(fe) does not become negligibly small below 
a given halo mass. However, in this case we can use the ap- 
proach introduced by Yoo et al. (2006): we write Xuik) as 
the sum of two terms, XM{k) = Xmi (k) + Xm^ {k), where: 

XAiAk) = f '""^ — Ui^{k\M) b{M) n{M) dM , 
Jo P 



Xm2 (k) 



M 

— Udm(fe|Af) b{M) n{M) dM . (24) 

A/„i„ P 



Following Yoo et al. (2006), we use the fact that «dm(fc|M) = 
1 over the relevant range of k as long as M is sufficiently 
small. This allows us to write 



Imi (fe) 



/* 

Jo 



M 



- b(M) n(M) dM 
P 

^ b(M) n(M) dM 
,„ P 



(25) 



where the last equality follows from the fact that the 
distribution of matter is by definition unbiased with re- 
spect to itself. Detailed tests have shown that this pro- 
cedure yields results that are sufficiently accurate as long 



Ha 



WMAP3 0.238 0.762 0.041 0.734 0.951 0.744 
WMAPl 0.3 0.7 0.040 0.7 1.0 0.9 



as Mmin < 10^° ft- 1 Me 
lO^ft-^M©. 



Throughout we adopt M„ 



2.3 Model Ingredients 

The computation of the galaxy-dark matter cross correlation 
(or its power spectrum) as outlined in the previous subsec- 
tions, requires the following ingredients: 

• The halo mass function, n{M), specifying the comoving 
number density of dark matter haloes of mass Al. 

• The halo bias function, fe(Af), which describes how 
haloes of mass M are biased with respect to the overall dark 
matter distribution. 

• The non-linear power spectrum of the dark matter dis- 
tribution, P<^(fc). 

• The mass density distribution of dark matter haloes, 
p{r\M). 

• The number density distribution of satellite galaxies in 
dark matter haloes, ns{r\M). 

• The halo occupation statistics for central and satellite 
galaxies, as parameterized by {Nc)m and {Ns}m, respec- 
tively. 

All these ingredients depend on cosmology. In this pa- 
per we consider two flat ACDM cosmologies. The first has a 
matter density Om = 0.238, a baryonic matter density Q^ ~ 
0.041, a Hubble parameter h = Ho/(100 kms"^ Mpc"^) = 
0.734, a power-law initial power spectrum with spectral in- 
dex n — 0.951 and a normalization ag = 0.744. These 
are the parameters that best-fit the 3-year data release of 
the Wilkinson Microwave Anisotropy Probe (WMAP, see 
Spergel et al. 2007), and we will refer to this set of cosmo- 
logical parameters as the WMAP3 cosmology. The second 
cosmology has Sim = 0.3, Q,h ~ 0.04, h — 0.7, n — 1.0 and 
(jg = 0.9. With strong support from the first year data re- 
lease of the WMAP mission (see Spergel et al. 2003), this 
choice of parameters has been considered in many previous 
studies. In what follows we will refer to this second set of 
parameters as the WMAPl cosmology. For clarity, the pa- 
rameters of both cosmologies are listed in Table 1. 

We define dark matter haloes as spheres with an average 
overdensity of 180, with a mass given by 



M=±L(180p)r?8o. 



(26) 



Here ri8o is the radius of the halo. We assume that dark 
matter haloes follow the NFW (Navarro, Frcnk & White 
1997) density distribution 



p(r) 



<5P 



{r/r,){l + r/r,y 



(27) 



where r« is a characteristic radius and 5 is a dimensionless 
amplitude which can be expressed in terms of the halo con- 
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centration parameter Cdi, 



; ri8o/r, as 



5 = 



180 4^ 

3 ln(l + Cdm) - Cdm/(1 + Cdm) 



(28) 



Numerical simulations show that Cdm is correlated with halo 
mass, and we use the relations given by Maccio et al. (2007), 
converted to our definition of halo mass. 

For the halo mass function, 'n{M), and halo bias func- 
tion, b{M), we use the functional forms suggested by War- 
ren et al. (2006) and Tinker et al. (2005), respectively, which 
have been shown to be in good agreement with numerical 
simulations. The linear power spectrum of density perturba- 
tions is computed using the transfer function of Eisenstein 
& Hu (1998), which properly accounts for the baryons, while 
the evolved, non-linear power spectrum of the dark matter, 
P^^{k), is computed using the fitting formula of Smith et 
al. (2003). 

For the number density distribution of the satellite 
galaxies, we assume a generalized NFW profile (e.g., van 
den Bosch et al. 2004): 



n.ir\M) « (^^ 



r 



(29) 



which scales as n^ oc r~° and Ub oc r~'^ at small and large 
radii, respectively. Similar to the dark matter mass distribu- 
tion, ns{r\M) has an effective scale radius TZr,, and can be 
parameterized via the concentration parameter Cg = Cdm/TZ. 
Observations of the number density distribution of satel- 
lite galaxies in clusters and groups seem to suggest that 
ns{r\M) is in good agreement with an NFW profile, for 
which a — 1 (e.g.. Beers & Tonry 1986; Carlberg, Yee & 
EUingson 1997a; van der Marel et al. 2000; Lin, Mohr & 
Stanford 2004; van den Bosch et al. 2005a). Several studies 
have suggested, however, that the satellite galaxies are less 
centrally concentrated than the dark matter, correspond- 
ing to 7^ > 1 (e.g., Yang et al. 2005; Chen 2007; More et 
al. 2008b). For our fiducial model we adopt a = TZ — 1, for 
which Us(A;|M) — Udin{k\M) (i.e., satellite galaxies follow 
the same number density distribution as the dark matter 
particles). In i^5.3l we examine how the results depend on a 
and 7?.. 

The final ingredient is a model for the halo occupation 
statistics. In their attempt to model the g-g lensing signal 
obtained from the SDSS, Seljak et al. (2005) and Mandel- 
baum et al. (2006) made the oversimplified assumption of 
a deterministic relation between central galaxy luminosity 
and host halo mass. In particular, they used 



{N,)MiLi,L2) 



1 if M = A'/(Li,L2) 
otherwise 



(30) 



where M{L\,L2) is the 'characteristic' mass of a halo that 
hosts a central galaxy with Li ^ L ^ L2. However, a realis- 
tic relation between central galaxy luminosity and host halo 
mass will have some scatter. As demonstrated by Tasitsiomi 
et al. (2004), this scatter can have an important impact on 
the g-g lensing signal (see also ^53] below). For the satellite 
galaxies, Seljak et al. (2005) and Mandelbaum et al. (2006) 
adopted a simple double power-law relation of the form 



{Ns}m(Li,L2) oc 



M 



ifAf ^ 3M(Li,L2) 
otherwise 



(31) 



In this paper we improve upon the analysis by Seljak et 
al. (2005) and Mandelbaum et al. (2006) by using a more re- 
alistic model for the halo occupation statistics. Furthermore, 
rather than fitting the model to the lensing data, we con- 
strain the occupation statistics using clustering data from 
the SDSS combined with a large galaxy group catalogue. 
Subsequently we use that model to predict the g-g lensing 
signal which we compare to g-g lensing data obtained from 
the SDSS. 

As a final remark, we emphasise that different quanti- 
ties, e.g. n{M), b{M), and PjJ^(fc), depend on redshift, z, 
even though we have not made this explicit in the equations. 



3 CONDITIONAL LUMINOSITY FUNCTION 

3.1 Model description 

In order to specify the halo occupation statistics, we use 
the CLF, <E>(I/|A/)dI/, which specifies the average number 
of galaxies with luminosities in the range L ± dL/2 that 
reside in a halo of mass M (Yang, Mo & van den Bosch 
2003; van den Bosch, Yang, Mo 2003). Following Cooray & 
Milosavljevic (2005) and Cooray (2006), we write the CLF 



<1>(LJM) = ML\M) + MMM) , 



(32) 



where <I>c(I/|Af) and "I>s(LjAf) represent central and satellite 
galaxies, respectively. The occupation numbers required for 
the computation of the galaxy-dark matter cross correlation 
then simply follow from 



{iVx)M(Ll,L2) 



<!>^{L\M)dL. 



(33) 



ii 



where 'x' refers to either 'c' (centrals) or 's' (satellites). Mo- 
tivated by the results of Yang, Mo & van den Bosch (2008; 
hereafter YMB08), who analyzed the CLF obtained from 
the SDSS galaxy group catalogue of Yang et al. (2007), we 
assume the contribution from the central galaxies to be a 
log-normal: 



-I'c(L|A/) 



27rln(10)crci 



exp 



(log L - log Lc 



2a2 



(34) 

Note that ctc is the scatter in log L (of central galaxies) at a 
fixed halo mass. Moreover, log Lc is, by definition, the expec- 
tation value for the (10-based) logarithm of the luminosity 
of the central galaxy, i.e.. 



log Lc = 



roo 
c = / <I>c 

Jo 



{L\M) logLdL. 



(35) 



For the contribution from the satellite galaxies we adopt a 
modified Schechter function: 



ML\M) 



exp 



(36) 



which decreases faster than a Schechter function at the 



bright end. Note that Lc, (Tc 



as and Li are all functions 



of the halo mass M. In the parametrization of these mass 
dependencies, we again are guided by the results of YMB08. 
In particular, for the luminosity of the central galaxies we 
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Table 2. Correlation lengths. 



Sample 


O-iM, - 


5 log ft 


< z> 


n) 


VI 


(-23.0, 


-21.5] 


0.173 


7.59 ± 0.75 


V2 


(-21.5, 


-21.0] 


0.135 


6.11 ±0.33 


V3 


(-21.0, 


-20.5] 


0.109 


5.62 ±0.16 


V4 


(-20.5, 


-20.0] 


0.089 


5.38 ±0.16 


V5 


(-20.0, 


-19.0] 


0.058 


4.90 ±0.18 


V6 


(-19.0, 


-18.0] 


0.038 


4.17 ±0.23 



Galaxy-galaxy clustering correlation lengths of Wang et al. (2007) 
used in this paper to constrain the CLF. Column (1) lists the ID 
of each volume limited sample, following the notation of Wang 
et al. (2007). Columns (2) and (3) indicate the absolute magni- 
tude range and the mean redshift of each sample, while column 
(4) lists the correlation length plus its standard deviation (in 
h~^ Mpc), obtained by fitting a power-law to the projected cor- 
relation function over the radial range [0.98, 9.6]/i~^ Mpc. See 
Wang et al. (2007) for details. 



adopt 



L4M) = Lo 



(M/Mi 



[1 + {M/Mi)Y' 



(37) 



so that Lc 
M > Ml. 



ex M^i for M < Ml and L^ ex M^^ fo^ 
Here Mi is a characteristic mass scale, and 
Lo = 2''^~'"'Lc{Mi) is a normalization. Using the SDSS 
galaxy group catalogue, YMB08 found that to good approx- 
imation 



L:(M) = 0.562ic(M) 



(38) 



and we adopt this parameterization throughout. For the 
faint-end slope and normalization of <I>s(LjM) we adopt 



as(M) = -2.0 + ai (1 arctan[o2 log(M/M2)] 1 , (39) 

and 

log[<^:(M)] = bo + &i(logMi2) + 62(logMi2)^ , (40) 

with Mi2 = M/(10^^/i"^ Mq). This adds a total of six free 
parameters: ai, 02, 6o, &i, &2 and the characteristic halo mass 
M2 . Neither of these functional forms has a physical motiva- 
tion; they merely were found to adequately describe the re- 
sults obtained by YMB08. Finally, for simplicity, and to limit 
the number of free parameters, we assume that (Tc(M) — CTc 
is a constant. As shown in More et al. (2008b), this assump- 
tion is supported by the kinematics of satellite galaxies in 
the SDSS. Thus, altogether the CLF has a total of eleven 
free parameters. 

Note that, with the parametrization of the CLF intro- 
duced above, the halo occupation statistics can be rewritten 
as: 



{N,)m{Li,L2) = I ' ML\M)dL = i[erf(a:2) 



erf(a;i 



(41) 



{N,)m{Li,L2) = 



2 ^2' 



/ $. 


{L\M)dL = 




Li 




LiV 
Li) 


-F 


2 +2' 


fL2 
Us 



(42) 



where erf(a::i) is the error function calculated at Xi = 
log(Z/i/I/c)/(\/2o"c) with i = 1,2 and F is the incomplete 
gamma function. 

As shown in Yang et al. (2003) and van den Bosch et 
al. (2003), the CLF can be constrained using the observed 
luminosity function, $(!/), and the galaxy-galaxy correla- 
tion lengths as a function of luminosity, ro{L). Here we 
use the luminosity function (hereafter LF) of Blanton et 
al. (2003a) uniformly sampled at 41 magnitudes covering the 
range -23.0 < "-^M^-Slog/i < -16.4. Here "'^Mr indicates 
the r-band magnitude K+E corrected to 2 = 0.1 following 
the procedure of Blanton et al. (2003b). For the correlation 
lengths as function of luminosity we use the results obtained 
by Wang et al. (2007) for six volume limited samples selected 
from the SDSS DR4. For completeness, these data are listed 
in Table 2. Finally, to strengthen our constraints, and to 
assure agreement with the CLF obtained from our SDSS 
group catalogue, we use the constraints on Lc{M), Qs(M) 
and (p'tiM) obtained by YMB08. 

For a given set of model parameters, we compute the 
LF using 



$(i) = / $(i|M) 
Jo 



n(M)dM. 



(43) 



The galaxy-galaxy correlation function for galaxies with lu- 
minosities in the interval [Li,L2] is computed using 



^es(r) =&gal(^l,^2)C(»")Cdm(r). 



(44) 



Here Cdm('') is the non- linear correlation function of the dark 



matter, which is the Fourier transform of P^itik), 



C(r) = 



[l + 1.17C^^(r)l^-^« 

[1+0.69C^L(,.)]2.09' 



(45) 



is the radial scale dependence of the bias as obtained by Tin- 
ker et al. (2005), and 6gai(-Li, L2) is the bias of the galaxies, 
which is related to the CLF according to 



bg!3.l{Ll, 1/2 ) = 



/p°°{iV)Mb(M)n(M)dM 
/(,°°(iV)Mn(M)dM 



(46) 



with 

{N)m = 



$(i|M)dL 



{N,)M{Li,L2) + {N,)MiLi,L2). 



(47) 



the average number of galaxies with luminosities in the range 
[Li, L2] that reside in a halo of mass M. 

To determine the likelihood function of our free pa- 
rameters we follow van den Bosch et al. (2007) and use 
the Monte-Carlo Markov Chain (hereafter MCMC) tech- 
nique. The goodness-of-fit of each model is judged using 



X 



xi + Xro + Xgc with 



xi 



2 

Xrg 



E 



$(LO-$(LO 



= E 



Cgg(7'0.0-l 

A4g(n),i) 



(48) 



(49) 
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Figure 2. Upper row, left and central panels. The luminosity function of galaxies and the luminosity dependence of the galaxy correlation 
length are plotted. Data come from the analysis of Blanton et al. (2003a) and Wang et al. 2007. The blue contours indicate the 68 and 95 
percent confidence level obtained from the MCMC. The agreement is extremely accurate for the luminosity function whereas is satisfying 
for the correlation length. Lower row, three panels. The additional information coming from the group catalogue of YMB08 is plotted 
together with the corresponding 68 and 95 percent confidence level derived with the MCMC. In particular, the halo mass dependence of 
the central galaxy luminosity, the satellite conditional luminosity function normalization 0* and the the exponent Og are shown in the 
left, central and middle panel, respectively. Upper row, right panel. The 68 and 95 percent confidence levels of the satellite fraction, /s, 
obtained from the CLF (see eq. 1511 ). 



Table 3. Best-fit CLF parameters obtained from SDSS clustering analysis 



Cosmology 
(1) 


logic 

(2) 


log Ml 
(3) 


71 
(4) 


72 
(5) 


ai 
(6) 


02 

(7) 


logM2 
(8) 


bo 
(9) 


bi 
(10) 


62 
(11) 


(12) 


(13) 


WMAP3 
WMAPl 


9.935 
9.915 


11.07 
11.16 


3.273 
3.336 


0.255 
0.248 


0.501 
0.484 


2.106 
2.888 


14.28 
14.54 


-0.766 
-0.854 


1.008 
0.906 


-0.094 
-0.062 


0.143 
0.140 


1.42 
1.70 



The best-fit CLF parameters obtained from the MCMC analysis for the WMAP3 and WMAPl cosmologies and the value of the 
corresponding reduced x^ ■ Masses and luminosities are in h~^ Mq and h~^ Lq, respectively. 



and 



2 
Xgc 



E 


9 

E 


9 

E 



\ogL,{Mi)-\ogL,{Mi) 
A log L, (MO 



,{Mi)-as{M^) 



A0*(MO 



(50) 



Here ? indicates an observed quantity and the subscripts 
'$', Vo' and 'GC refer to the luminosity function, the 
galaxy-galaxy correlation length and the group catalogue, 



respectively. Note that, by definition, ^gg(ro,i) = 1. Table 4 
lists the best-fit parameters obtained with the MCMC tech- 
nique for both the WMAPl and WMAP3 cosmologies, as 
well as the corresponding value of Xred = X^ /Ndoi. Here 
Ndoi = 74 — 11 = 63 is the number of degrees of freedom. 



3.2 Results 

Fig.[2]shows the results obtained for the WMAP3 cosmology. 
In each panel the blue contours indicate the 68 and 95 per- 
cent confidence levels obtained from the MCMC. The upper 
left-hand panel shows that the CLF model accurately fits 
the galaxy LF of Blanton et al. (2003a). The fit to the cor- 
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Figure 3. Same as Fig.[2]but for the WMAPl cosmology. 



relation lengths as function of luminosity, shown in the up- 
per middle panel, is less accurate, although data and model 
typically agree at the Ict level. The lower panels of Fig. [2] 
show the 68 and 95 percent confidence levels on Lc{M), 
cj>t{M) and Qs(Af), compared with the results obtained by 
YMB08 from the SDSS group catalogue of Y07. Since these 
data have been used as additional constraints, it should not 
come as a big surprise that the CLF is in good agreement 
with these data. We emphasise, though, that it is not triv- 
ial that a single halo occupation model can be found that 
can simultaneously fit the LP, the luminosity dependence 
of the galaxy-galaxy correlation functions, and the results 
obtained from a galaxy group catalogue. 

Finally, the upper right-hand panel of Fig. [2] shows the 
satellite fraction. 



UL) = 



$(L) 



$s(L|M)n(M)dM. 



(51) 



as function of luminosity. This is found to decrease from 
~ 0.27 ± 0.03 at °-^Mr - 51og/i = -17 to virtually zero at 
^'^Mr — 51og/i = —23. The fact that the satellite fraction 
decreases with increasing luminosity is in qualitative agree- 
ment with previous studies (Tinker et al. 2006; Mandelbaum 
et al. 2006; van den Bosch et al. 2007). 

We have repeated the same exercise for the WMAPl 
cosmology. As evident from Pig. O for this cosmology we 
can obtain a CLP that fits the data almost equally well (the 
reduced x^ is only slightly higher than for the WMAP3 cos- 
mology; see Table 3). Note that the group data (shown in 
the lower panels) differ from that in Pig. [2j even though the 
group catalogue is the same. This owes to the fact that the 



halo mass assignments of the groups are cosmology depen- 
dent (see Y07 for details). The satellite fractions inferred 
for this cosmology, shown in the upper right-hand panel 
of Pig. [31 are similar, though slightly higher, than for the 
WMAP3 cosmology, in excellent agreement with van den 
Bosch et al. (2007). 

The fact that both cosmologies allow an (almost) 
equally good fit to these data, despite the relatively large 
differences in halo mass function and halo bias, illustrates 
that the abundance and clustering properties of galaxies al- 
low a fair amount of freedom in cosmological parameters. 
However, as demonstrated in van den Bosch, Mo & Yang 
(2003), the best-fit CLPs for different cosmologies predict 
different mass-to-light ratios as function of halo mass. This 
is evident from Pig. U which shows the mass-to-light ratios 
Af/{Li9,5) A/ as function of halo mass inferred from our CLP 
MCMCs for the WMAPl and WMAP3 cosmologies. Here 
{^19.5)a/ is the average, total luminosity of all galaxies with 
"■^Mr — 5 log ft ^ —19.5 that reside in a halo of mass M, 
which follows from the CLP according to 



(il9.5)j\/ = 



$(L|A/)LdL, 



(52) 



with Lmi-n the luminosity corresponding to a magnitude 
"■^Afr — 51og/i — —19.5. Clearly, the mass-to-light ratios 
inferred for the WMAPl cosmology are significantly higher 
than for the WMAP3 cosmology (see also van den Bosch 
et al. 2007, where a similar result was obtained using data 
from the 2dPGRS). Hence, the abundance and clustering 
properties of galaxies can be used to constrain cosmological 
parameters, as long as one has independent constraints on 
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Figure 4. The 68 and 95 percent confidence levels for the mass- 
to-light ratios, M/{Li9.5>m, obtained from the CLF MCMCs for 
the WMAP3 and WMAPl cosmologies. 



the mass-to-light ratios of dark matter haloes. This is ex- 
actly what is provided by g-g lensing. In the next section, 
we therefore use the CLF models presented here to predict 
the g-g lensing signal, which we compare to SDSS data. 



4 GALAXY-GALAXY LENSING 
4.1 Model Predictions 

In order to compute the BSD, AE, as a function of the co- 
moving separation on the sky, R, we proceed as follows. We 
start by calculating the four different terms of the galaxy- 
dark matter cross power spectrum defined in i; l2.1l and ij 12.21 
Next we inverse Fourier transform each of these terms using 



C'Ur) 



1 

2^ 



/.P,%(fc)^fcMfc, (53) 



where '/i' stands for Ih or 2h, and 'x' refers to either 'c' 
(centrals) or 's' (satellites). These are used to compute the 
corresponding four terms of the surface density, E^'^(_R), 



{R) = 2-p ed^(r) 



rdr 



Vr^ 



i?2 



(54) 



Note that we are allowed to neglect the contribution coming 
from the constant background density, p, (cf. equation [5]) 
because it will cancel in defining the ESD (this is known in 
gravitational lensing theory as the mass-sheet degeneracy). 
The final ESD then simply follows from 



AE(i?) = AE''^'"(i?) + AE"^'''(i?) 
+ AE"^'"(i?)-FAE"^''=(i?). 



(55) 



in which the relative weight of each term is already included 
via the central and satellite fractions in the definitions of 
the corresponding power spectrslj. 

Before comparing the g-g lensing predictions from our 
CLF models to actual data, we first demonstrate how the 
four different terms contribute to the total ESD. The left- 
hand panel of Fig. [S] shows the AE(7?) obtained from our 
best-fit CLF model for the WMAP3 cosmology for three dif- 
ferent luminosity bins, as indicateqfl. Note that the fainter 
luminosity bins reveal a more 'structured' excess surface 
density profile, with a pronouced 'bump' at _R ~ lh~^ Mpc, 
which is absent in the AE(_R) of the brighter galaxies. The 
reason for this is evident from the middle and right-hand 
panels of Fig. [5] which show the contributions to AE(_R) 
from the four different terms for the faint ( — 16 ^ "'^Af,. — 
51og/i ^ -17) and bright (-21 > °-^Af^ - 51og/i > -22) 
luminosity bins, respectively. In both cases, the 1-halo cen- 
tral term dominates on small scales. In the case of the faint 
galaxies, the 1-halo satellite term dominates over the radial 
range 0.1/i~^ Mpc <, _R <, 5/i~^ Mpc, and is responsible for 
the pronounced bump on intermediate scales. In the case 
of the bright bin, however, the 1-halo central term domi- 
nates all the way out to 7? ~ 2h~^ Mpc. This owes to the 
fact that bright centrals reside in more massive haloes, which 
are larger and cause a stronger lensing signal, and due to the 
fact that the satellite fraction of brighter galaxies is smaller. 
The fact that the 1-halo satellite term peaks at intermediate 
scales, rather than at i? = 0, owes to the fact that AE^'^'°(7?) 
reflects a convolution of the host halo mass density pro- 
flle with the number density distribution of satellite galax- 
ies. On large scales {R >, 3/i~^ Mpc), which roughly reflects 
two times the virial radius of the most massive dark mat- 
ter haloes, the ESD is dominated by the 2-halo terms. Note 
that the faint galaxies, with —16 ^ ^'^Mr — 51og/i ^ —17, 
have the same large scale ESD as the bright galaxies with 
—21 ^ "'^Mr- 51og/i ^ —22, indicating that they have sim- 
ilar values for their bias. This owes to the fact that many 
of the faint galaxies are satellites which reside in massive 
haloes. Note also that the 2-halo central term reveals a 
fairly abrupt truncation at small radii, which owes to halo- 
exclusion (see Appendix). This truncation also leaves a sig- 
nature in the total lensing signal, which is more pronounced 
for the fainter lenses. We caution, however, that the sharp- 
ness of this feature is partially an artefact due to our ap- 
proximate implementation of halo-exclusion. Nevertheless, 
it is clear from Fig. [5] that the excess surface densities ob- 
tained from g-g lensing measurements contain a wealth of 
information regarding the galaxy-dark matter connection. 



4.2 Data 

The g-g lensing data used here is described in Seljak et 
al. (2005) and Mandelbaum et al. (2006) and has been 
kindly provided to us by R. Mandelbaum. The catalogue of 



^ Note that our notation differs slightly from that in Mandelbaum 
et al. (2006) 

^ Here, for simplicity, we have used the halo mass function and 
halo bias function computed at 2 = 0. In i|4.2l when we compare 
our models to data, we will use the halo mass function and halo 
bias function at the average redshift of the lenses instead. 
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Figure 5. The predicted ESD up to large scales (i? ~ 30/i~^Mpc) for three luminosity bins, as indicated. The solid lines refer to the 
total signal as predicted according to our model. The dotted lines refer to the 1-halo central term, whereas the dashed lines refer to the 
1-halo satellite term. Note that they dominate the signal on different scales (see text). The long dashed lines refer to the 2-halo central 
term. It rises steeply at relatively large scales due to our halo exclusion treatment (see Appendix). The 2-halo satellite term is indicated 
with the dotted-dashed line. 



Table 4. Luminosity bins of the SDSS g-g lensing data 



ID 


■^Mr-^\ogh 


(^> 


(1) 


(2) 


(3) 


LI 


-18.0,-17.0] 


0.032 


L2 


-19.0,-18.0] 


0.047 


L3 


-20.0, -19.0] 


0.071 


L4 


-21.0,-20.0] 


0.10 


L5f 


-21.5,-21.0] 


0.14 


L5b 


-22.0,-21.5] 


0.17 


L6f 


-22.5,-22.0] 


0.20 



Luminosity bins of the lenses. Column (1) lists the ID of each lu- 
minosity bin, following the notation of Mandelbaum et al. (2006). 
Column (2) indicates the magnitude range of each luminosity bin 
(all magnitudes are K+E corrected to z = 0.1). Column (3) indi- 
cates the mean redshift of the lenses in each luminosity bin. See 
Mandelbaum et al. (2006) for details. 

lenses consists of 351, 507 galaxies with magnitudes —17 ^ 
°-^Mr - 51og/i > -23 and redshifts 0.02 < z < 0.35 taken 
from the main galaxy catalogue of the SDSS Data Release 
4 (Adelman-McCarthy et al. 2006). This sample is split in 
7 luminosity bins (see Table 4), and for each of these lumi- 
nosity bins the excess surface density profiles, AE(7?), have 
been determined from measurements of the shapes of more 
than 30 million galaxies in the SDSS imaging data down to 
an apparent r-band magnitude of r = 21.8. The resulting 
data are shown as solid dots with errorbars in Fig. [S] We 
refer the reader to Mandelbaum et al. (2005b, 2006) for a 
detailed description of the data and of the methods used to 
determine the ESD profiles. 



4.3 Results for the WMAP3 Cosmology 



Using the methodology outlined in g2] and il4.1l we now use 
the CLF for the WMAP3 cosmology obtained in ^ to pre- 



dict the g-g lensing signal for the 7 luminosity bins listed in 
Table 4. For each luminosity bin we compute the ESD pro- 
file, AE(iZ), at the mean redshift of the sample, i.e., we use 
the halo mass function, n[M), the halo bias function, b{M), 
and the non- linear power spectrum, P^rtik) that correspond 
to the mean redshift listed in the third column of Table 4. 
We have verified, though, that computing AE(i?) simply at 
z = instead has a negligible impact on the results. 

The results are shown in Fig. |S1 where the solid dots 
with errorbars correspond to the SDSS data and the solid 
lines are the predictions of our best-fit CLF model (whose 
parameters are listed in Table 3). Note that this model fits 
the data remarkably well, which is quantified by the fact 
that the reduced x^ is 3.1. We emphasize that there are no 
free parameters here: the ESD has been computed using the 
CLF that has been constrained using the LF and the clus- 
tering data. The good agreement between model and lensing 
data thus provides independent support for the halo occu- 
pation statistics described by our WMAP3 CLF model, in 
particular for the mass-to-light ratios and satellite fractions, 
which have an important impact on the lensing signal. 

The different curves in each of the panels in Fig. [6] 
show the contribution to the lensing signal due to the four 
separate terms: AE^'^''' (dotted lines), AE^''''' (short-dashed 
lines), AE^'"''^ (long-dashed lines), and AE^'''^ (dot-dashed 
lines). In agreement with the examples shown in Fig. \S\ the 
1-halo central term becomes increasingly more dominant for 
more luminous lenses. In fact, in the brightest luminosity 
bin (L6f) it dominates over the entire radial range probed. 
In the low-luminosity bins, most of the observed lensing sig- 
nal at i? >, 200h~^ kpc is dominated by the 1-halo satellite 
term. The fact that our model accurately fits the data, thus 
supports the satellite fractions inferred from our CLF model, 
and shown in the upper right-hand panel of Fig. (2] 

Both Seljak et al. (2005) and Mandelbaum et al. (2006) 
did not account for the contributions of the 2-halo terms in 
their analyses of the galaxy-galaxy lensing signal. Our model 
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Figure 6. The excess surface density AS as a function of the comoving transverse separation R is plotted for different bins in luminosity 
of the lens galaxy (see Table 4). The solid line represents the total signal as predicted by the model, data points and error bars come from 
Seljak et al. (2005), see text. The different contributions to the signal are also plotted. The dotted line represents the 1-halo central term 
which obviously dominates at the smallest scales in all cases. Note that this term dominates on larger and larger scales when brighter 
galaxies are considered, reflecting the idea that brighter galaxies live on average in more massive haloes. The dashed line represents 
the 1-halo satellite term which is dominant only for faint galaxies and only on intermediate scales (0.1-1 h~^Mpc). The 2-halo central 
is plotted with a long dashed line and it becomes relevant on large scales {R > l/i~^Mpc). Note that the strong truncation of this 
term at small scales is due to our implementation of halo exclusion (see Appendix). The 2-halo satellite term (dotted-dashed line) never 
dominates but it can contribute up to 20% of the total signal. 



indicates that, although the 2-halo terms never dominate 
the total signal, they can contribute as much as 50 percent 
at large radii {R ~ lh~^ Mpc). We thus conclude that the 
2-halo terms cannot simply be ignored. 



4.4 Comparison with the WMAPl Cosmology 



As shown in 93^ the WMAP3 and WMAPl cosmologies 
both allow a good fit to the clustering data, luminosity func- 
tion and galaxy group results. However, the corresponding 
CLFs predict mass-to-light ratios that are significantly dif- 
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Figure 7. The predictions for the lensing signal, AS(i?), are shown for two different sets of cosmological parameters (WMAPl and 
WMAP3, see text). The green (blue) shaded area corresponds to the 95% confidence level of the WMAPl (WMAP3) model. Note that, 
although the cosmological parameters of these two cosmologies only differ by < 20 percent (see Table 1), the ESD predictions are very 
different, and can easily be discriminated. 



ferent. Since the galaxy-galaxy lensing signal is very sensi- 
tive to these mass-to-light ratios, it is to be expected that 
our WMAP3 and WMAPl CLFs will predict significantly 
different ESD profiles, thus allowing us to discriminate be- 
tween these two cosmologies. 

Fig.[7]shows the 95 percent confidence levels for AE(_R) 
obtained from our CLF MCMCs for both the WMAP3 
(blue) and WMAPl (green) cosmologies. Indeed, as antici- 
pated, for the WMAPl cosmology we obtain excess surface 
densities that are significantly higher than for the WMAP3 



cosmology, in accord with the higher mass-to-light ratios (cf. 
Fig. ID). A comparison with the SDSS data clearly favors the 
WMAP3 cosmology over the WMAPl cosmology. In fact, for 
the latter our best-fit CLF model yields a reduced x^ of 29.5, 
much larger than for the WMAP3 cosmology (x?od — 3.1). 
Note that the cosmological parameters for these two cos- 
mologies are very similar: flm and ag differ only by ~ 20 
percent (in addition to a ~ 5 percent difference in n). Yet, 
we can very significantly favor one cosmology over the other. 
This indicates that the combination of clustering data and 
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Figure 8. The average number of central galaxies as a function 
of halo mass obtained from our best-fit CLF for the WMAP3 
cosmology. This is equivalent to the probability VciMlLi, L2) 
that a central galaxy with Li ^ L ^ L2 is hosted by a halo 
of mass M. Results are shown for four different luminosity bins, 
as indicated. Note that brighter centrals reside, on average, in 
more massive haloes. In addition, the width of 'Pc(M\L\, L2) also 
increases with luminosity. 



g-g lensing data can be used to put tight constraints on cos- 
mological parameters. A detailed analysis along these lines 
is deferred to a forthcoming paper (Cacciato et al. , in prepa- 
ration). 



5 MODEL DEPENDENCIES 

Although our computation of the g-g lensing signal does not 
involve any free parameters (these are already constrained 
by the clustering data), a number of assumptions are made. 
In particular, haloes are assumed to be spherical and to 
follow a NFW density distribution, central galaxies are as- 
sumed to reside exactly at the center of their dark matter 
haloes, and satellite galaxies are assumed to follow a ra- 
dial number density distribution that has the same shape 
as the dark matter mass distribution. In addition, we made 
assumptions regarding the functional form of the CLF. Al- 
though most of these simplifications are reasonable, and 
have support from independent studies, they may have a 
non-negligible impact on the predictions of the g-g lensing 
signal. If this is the case, they will affect the reliability of 
the cosmological constraints inferred from the data. In this 
section we therefore investigate how strongly our model pre- 
dictions depend on these oversimplified assumptions. 

Some of these dependencies were already investigated 
in our companion paper (Li et al. 2008). In particular, we 
have shown that the fact that realistic dark matter haloes 
are ellipsoidal, rather than spherical, can be safely ignored 
(i.e., its impact on the ESD profiles is completely negligible). 
On the other hand, if central galaxies are not located exactly 



at the center of their dark matter haloes, this may have a 
non-negligible impact on the lensing signal on small scales 
{R <, 0.1h~^ Mpc). Fortunately, for realistic amplitudes of 
this offset (van den Bosch et al. 2005b), the effect is fairly 
small and only restricted to the most luminous galaxies (see 
Li et al. 2008 for details). 

Below we investigate three additional model dependen- 
cies: the scatter in the relation between hght and mass, the 
concentration of dark matter haloes, and the radial number 
density distribution of satellite galaxies. To that extent we 
compare our fiducial model (the best-fit CLF model for the 
WMAP3 cosmology presented above), to models in which 
we change only one parameter. 

5.1 Scatter in the Lc — M relation 

An important improvement of our halo occupation model 
over that used by Seljak et al. (2005) and Mandelbaum et 
al. (2006) is that we allow for scatter in the relation be- 
tween hght and mass. In particular, we model the proba- 
bility function 'Pc{L\M) = <^c{L\M) as a log-normal with a 
scatter, CTc, that is assumed to be independent of halo mass. 
As demonstrated in More et al. (2008b), this assumption is 
consistent with the kinematics of satellite galaxies, and it is 
supported by semi-analytical models for galaxy formation. 
Note, though, that this does not imply that the scatter in 
Pc(M|L), which is the probability function which actually 
enters in the computation of the g-g lensing signal, is also 
constant. In fact, it is not. This is illustrated in Fig.[8l which 
shows 'Pc{M\L\,L2) of our fiducial model for four luminos- 
ity bins. Two trends are evident: more luminous centrals 
reside, on average, in more massive haloes, and they have 
a larger scatter in halo masses. As discussed in More, van 
den Bosch & Cacciato (2008), the fact that the scatter in 
Vc[M\L) increases with luminosity simply owes to the fact 
that the slope of Lc{M) becomes shallower with increasing 
M (see the lower right-hand panels of Figs. [2] and [3} . As is 
evident from Fig. [HI this is a strong effect, with the scatter in 
■Pc(M|L) becoming extremely large at the bright end. Note 
that this scatter is not dominated by the width of the lumi- 
nosity bin. Hence, even if one were able to use infinitesimally 
narrow luminosity bins when stacking lenses, the scatter in 
Pc(Af jL) will still be very appreciable. 

As first shown by Tasitsiomi et al. (2004) , scatter in the 
relation between light and mass can have a very significant 
impact on the ESDs. This is demonstrated in the upper pan- 
els of Fig. [9] which show the impact on AE(7?) of changing 
<Jc by 0.05 compared to our best-fit CLF value of ctc = 0.14; 
all other parameters are kept fixed at their fiducial values 
(see Table 3). Note that these changes in CTc have a negli- 
gible impact on AE(_R) for the low luminosity bins. At the 
bright end, however, relatively small changes in ctc have a 
very significant impact on AE(_R). In particular, increas- 
ing the amount of scatter reduces the ESD. This behavior 
owes to the shape of the halo mass function. Increasing the 
scatter adds both low mass and high mass haloes to the dis- 
tribution, and the overall change in the average halo mass 
depends on the slope of the halo mass function. Brighter 
galaxies live on average in more massive haloes where the 
halo mass function is steeper. In particular, when the aver- 
age halo mass is located at the exponential tail of the halo 
mass function, an increase in the scatter adds many more 
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Figure 9. The impact of various model parameters on AE(H). Results are shown for three luminosity bins, as indicated at the top of each 
column. In each panel the solid line corresponds to our fiducial model (the best-fit CLF model for the WMAP3 cosmology presented in 
Fig-O, while the dotted and dashed lines correspond to models in which we have only changed one parameter or model ingredient. Upper 
panels: the impact of changes in the parameter ctc, which describes the amount of scatter in <3?c(-L|M) (see equation [34]). Second row 
from the top: the impact of changes in the halo concentration, c^^{M). In particular, we compare three models for the mass dependence 
of Cdm; Maccio et al. (2007; MAC), Bullock et al. (2001; BUL), and Eke et al. (2001; ENS). Third row from the top: the impact of 
changes in 7?. = c^jm /cs , which controls the concentration of the radial number density distribution of satellite galaxies relative to that of 
the dark matter. Lower panels: The impact of changes in a, which specifies the central slope of the radial number density distribution 
of satellite galaxies. See text for a detailed discussion. 



low mass haloes than massive haloes, causing a drastic shift 
in the average halo mass towards lower values. On the other 
hand, fainter galaxies live in less massive haloes, where the 
slope of the halo mass function is much shallower. Conse- 



quently, a change in the scatter does not cause a significant 
change in the average mass. 

Clearly, if the g-g lensing signal is used to constrain cos- 
mological parameters, it is important that one has accurate 
constraints on (Jc. From the clustering analysis presented in 
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^3:2] we obtain 0.14±0.01 (for both WMAPl and WMAP3). 
This is in good agreement with previous studies: Cooray 
(2006), using a CLF to model the SDSS r-band LF, obtained 



0.171 



YMB08, using a SDSS galaxy group cata- 



logue, obtained a^ = 0.13 ± 0.03, and More et al. (2008b), 
using the kinematics of satellite galaxies in the SDSS, find 
(Jc = 0.16 ± 0.04 (all errors are 68% confidence levels). Al- 
though it is reassuring that very different methods obtain 
values that are in such good agreement, it is clear that the 
remaining uncertainty may have a weak impact on our abil- 
ity to constrain cosmological parameters. Fortunately, the 
scatter only impacts the results at the bright end, so that 
one can always check the results by removing data from the 
brightest luminosity bins. 



5.2 The dark matter halo concentration 

The g-g lensing signal on small scales reflects the projected 
mass distribution of the haloes hosting the lensing galaxies. 
Therefore, the detailed shape of AE(_R) on small scales is 
sensitive to the mass distribution of dark matter haloes. In 
our model, we have assumed that dark matter haloes follow 
NFW profiles, which are characterized by their concentra- 
tion parameters, Cdm- Halo concentrations are known to de- 
pend on both halo mass and cosmology, and various analyti- 
cal models have been developed to describe these dependen- 
cies (Navarro et al. 1997; Bullock et al. 2001; Eke, Navarro & 
Steinmetz 2001; Neto et al. 2007; Maccio et al. 2007, 2008). 
Unfortunately, these models make slightly different predic- 
tions for the mass dependence of Cdm (mainly due to the fact 
that the numerical simulations used to calibrate the models 
covered different limited mass ranges). In Li et al. (2008), 
we have shown that changing Cdm by a factor of two has a 
very large impact on the BSD profiles. However, this is much 
larger than the typical discrepancies between the different 
models for Cdm(M). The second row of panels in Fig.[9]shows 
AE(_R) obtained for three of these models: the solid lines (la- 
belled MAC) corresponds to our fiducial model for which we 
have used the Cdm(A^) relation of Maccio et al. (2007). The 
dotted lines (labelled BUL) and dashed lines (labelled ENS) 
correspond to the Cdm(Af) relations of Bullock et al. (2001) 
and Eke et al. (2001), respectively. The BUL model predicts 
halo concentrations that are about 15 percent higher than 
for the MAC model. The ENS model predicts a c^^{M) that 
is somewhat shallower than the BUL and MAC models. As 
is evident from Fig. [9] though, the results based on these 
three different models are very similar. We thus conclude 
that our results are robust to uncertainties in the relation 
between halo mass and halo concentration. 



5.3 Number density of satellite galaxies 

In our modelling of the g-g lensing signal, we have assumed 
that the number density distribution of satellite galaxies can 
be described by a generalised NFW profile (eq. [29]), which 
is parameterized by two free parameters: a and TZ. In the 
models presented above, we have assumed that a = TZ = 1, 
so that the number density distribution of satellite galax- 
ies has exactly the same shape as the dark matter distri- 
bution. As discussed in ii2.3l though, there is observational 
evidence which suggests that satellite galaxies are spatially 



anti- biased with respect to the dark matter (i.e., their radial 
distribution is less concentrated than that of the dark mat- 
ter). This is also supported by numerical simulations, which 
show that dark matter subhaloes (which are believed to host 
satellite galaxies) are also spatially anti-biased with respect 
to the dark matter (e.g., Moore et al. 1999, De Lucia 2004). 

The panels in the third row of Fig. [5] show the impact 
of changing the concentration of the radial number density 
distribution of satellite galaxies. In particular, we compare 
the ESD profiles obtained for our fiducial model {TZ = 1.0, 
solid lines) with models in which TZ = 0.5 (dotted lines) and 
TZ = 2.Q (dashed lines). Recall that TZ = Cdm/cs, so that 
TZ > 1 {TZ < 1) corresponds to satellite galaxies being less 
(more) centrally concentrated than the dark matter. Note 
that changes in 7?. have a negligible effect on AE(i?) for the 
bright luminosity bins. This simply owes to the fact that the 
ESD of bright lenses is completely dominated by the 1-halo 
central term (i.e., the satellite fraction of bright galaxies is 
very small). For the fainter luminosity bins, however, an in- 
crease (decrease) in 7i causes a decrease (increase) in AS(7?) 
on intermediate scales (0.1/i~^Mpc ^ R ^ lh~^Mpc), which 
is the scale on which the 1-halo satellite term dominates. The 
effect, though, is fairly small (typically smaller than the er- 
rorbars on the data points). 

The last row of Fig. [5]shows the impact of changing the 
central slope, q, of ns{r). If the number density distribu- 
tion of satellite galaxies has a central core (a = 0), rather 
than a NFW- like cusp {a — 1), it has a similar impact on 
the lensing signal as assuming a less centrally concentrated 
ns{r). In fact, the ESD profiles for (a,7^) = (0.0,1.0) are 
very similar to those for {a,TZ) = (1.0, 2.0). The main con- 
clusion, though, is that our results are not very sensitive to 
the exact form of Us (r) (see also Yoo et al. 2006). Clearly, our 
conclusion that the WMAP3 cosmology is strongly preferred 
over the WMAPl cosmology is not affected by uncertainties 
in the radial distribution of satellite galaxies. 



6 CONCLUSIONS 

Galaxy clustering and galaxy-galaxy lensing probe the 
galaxy-dark matter connection in complementary ways. 
Since the clustering of dark matter haloes depends on cos- 
mology, the halo occupation statistics inferred from the ob- 
served clustering properties of galaxies are degenerate with 
the adopted cosmology. Consequently, different cosmologies 
imply different mass-to-light ratios for dark matter haloes. 
Galaxy-galaxy lensing, on the other hand, yields direct con- 
straints on the actual mass-to-light ratios of dark matter 
haloes. Combined, clustering and lensing therefore offer the 
opportunity to constrain cosmological parameters. 

Although the advent of wide and deep surveys has re- 
sulted in clear detections of galaxy-galaxy lensing, a proper 
interpretation of these data in terms of the link between 
galaxies and dark matter haloes has been hampered by the 
fact that the lensing signal can only be detected when stack- 
ing the signal of many lenses. Since not all lenses reside 
in haloes of the same mass, the resulting signal is a non- 
trivial average of the lensing signal due to haloes of different 
masses. In addition, central galaxies (those residing at the 
center of a dark matter halo) and satellite galaxies (those 
orbiting around a central galaxy) contribute very different 
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lensing signals, even when they reside in haloes of the same 
mass (e.g., Yang et al. 2006). This has to be properly ac- 
counted for, and requires knowledge of both the satellite 
fractions and of the spatial number density distribution of 
satellite galaxies within their dark matter haloes. 

In this paper, we model galaxy-galaxy lensing with the 
CLF, $(I/|M), which describes the average number of galax- 
ies of luminosity L that reside in a halo of mass M. This CLF 
is ideally suited to model galaxy-galaxy lensing. In particu- 
lar, it is straightforward to account for the fact that there 
is scatter in the relation between the luminosity of a central 
galaxy and the mass of its dark matter halo. This repre- 
sents an improvement with respect to previous attempts to 
model the g-g lensing signal obtained from the SDSS, which 
typically ignored this scatter (e.g. Seljak et al. 2005; Man- 
delbaum et al. 2006). However, in agreement with Tasitsiomi 
et al. (2004), we have demonstrated that the scatter in this 
relation has an important impact on the g-g lensing signal 
and cannot be ignored. We also improved upon previous 
studies by modelling the 2-halo term (the contribution to 
the lensing signal due to the mass distribution outside of 
the halo hosting the lens galaxy), including an approximate 
treatment for halo exclusion. 

Following Cooray & Milosavljevic (2005), we split the 
CLF in two components: one for the central galaxies and 
one for the satellites. This facilitates a proper treatment of 
their respective contributions to the g-g lensing signal. The 
functional forms for the two CLF components are motivated 
by results obtained by Yang et al. (2008) from a large galaxy 
group catalogue. For a given cosmology, the free parameters 
of the CLF are constrained using the luminosity function, 
the correlation lengths as function of luminosity, and some 
properties extracted from the group catalogue. We have per- 
formed our analysis for two different ACDM cosmologies: the 
WMAPl cosmology, which has fim = 0.3 and erg = 0.9 and 
the WMAP3 cosmology with fi^ = 0.238 and crs = 0.744. 
For both cosmologies we have obtained CLFs that can accu- 
rately fit the abundances and clustering properties of SDSS 
galaxies. However, these CLFs predict mass-to-light ratios 
that are very different. This reflects the degeneracy between 
cosmology and halo occupation statistics alluded to above. 
In order to break this degeneracy, we use these CLFs to 
predict the g-g lensing signal (with no additional free pa- 
rameters), which is compared to the SDSS data obtained by 
Seljak et al. (2005) and Mandelbaum et al. (2006). While the 
WMAP3 CLF predictions are in excellent agreement with 
the data, the CLF for the WMAPl cosmology predicts ex- 
cess surface densities that are much higher than observed. 
Although the cosmological parameters of the WMAPl and 
WMAP3 cosmologies only differ at the 20 percent level, the 
combination of clustering and lensing allows us to strongly 
favor the WMAP3 cosmology over the WMAPl cosmology. 
In a companion paper by Li et al. (2008), we use a com- 
pletely different technique to model g-g lensing, but never- 
theless reach exactly the same conclusion. 

In order to test the robustness of our results we have 
performed a number of tests. In particular, we have shown 
that small uncertainties in the expected concentrations of 
dark matter haloes, or in the radial number density distri- 
butions of satellite galaxies, only have a very small impact on 
the predicted lensing signal. In addition, although our treat- 
ment of halo exclusion is only approximate, we have demon- 



strated that it is sufficiently accurate. Finally, as shown by 
Li et al. (2008), making the oversimplified assumption that 
dark matter haloes are spherical rather than ellipsoidal also 
has a negligible impact on the lensing predictions. We thus 
conclude that our method yields accurate and reliable pre- 
dictions for g-g lensing. 

To summarize, as already discussed by Yoo et al. (2006) , 
the combination of clustering and lensing can be used to 
put tight constraints on cosmological parameters. In this pi- 
lot study we have shown that current data from the SDSS 
strongly favors the WMAP3 cosmology over the WMAPl 
cosmology. In a follow-up paper we will present a more de- 
tailed analysis of the cosmological constraints that can be 
obtained using this technique. 
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APPENDIX A: HALO EXCLUSION 

By definition, the 2-halo terms of the galaxy-dark matter 
cross correlation, ^g,dm(?"), only considers pairs of galaxies 
and dark matter particles that reside in different haloes. 
Since two haloes can not overlap spatially, this implies that 
the 2-halo terms given by Eqs. (|20|l and (|21[) need to be 
modified to take account of halo exclusion. The concept of 
halo exclusion is illustrated in Eig. I All for the 2-halo central 
and 2-halo satellite terms separately. Consider a spherical 
halo of mass M and radius ri8o that hosts a central galaxy. 
It is clear that this central galaxy cannot contribute any 
signal to the 2-halo term on scales smaller than ri8o. Hence, 
if all central galaxies lived in haloes of a fixed mass M, 
then 1 + Cg dm('") = for r < ri8o. In reality, though, one 
needs to account for the fact that centrals occupy haloes 
that span a range in halo masses, even if the centrals all 
have the same luminosity. In the case of the satellite galaxies 
the situation is even more complicated. In particular, since 
satellite galaxies can reside at the outskirts of dark matter 
haloes, the 2-halo satellite term can still have non-zero power 
at r <C risQ. Thus, halo exclusion has less impact on the 2- 
halo satellite term than on the 2-halo central term. 

Although the concept of halo exclusion is quite simple, 
a proper implementation of it in the halo model is extremely 
tedious numerically. We therefore use only an approximate 
treatment, which has the advantage that it is straightfor- 
ward to implement numerically. Eirst of all, we ignore halo 
exclusion for the 2-halo satellite term. Since this term is al- 
ways smaller than the 2-halo central term, and since halo 
exclusion is less important for satellites than for centrals, 
this should not have a significant impact on the results. Eor 
the 2-halo central term we proceed as follows: for each lu- 
minosity bin, [Li,L2], we simply set 1 -I- ^g_dm(^) ~ f^'' 
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Figure Al. Illustration of halo exclusion. The upper panel shows 
two haloes, of masses M and A/', and corresponding radii rigo 
and J'igo' respectively. The halo of mass M hosts a central galaxy. 
Since two haloes cannot overlap, this central galaxy does not con- 
tribute any signal to the 2-halo central term of the galaxy-dark 
matter cross correlation function on scales r < rigo- In the case 
of the 2-halo satellite term, illustrated in the lower panel, there is 
still a contribution even on very small scales (r <C ^180 )i simply 
because satellite galaxies can reside near the edge of the halo. 
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on small scales is completely dominated by the l-halo terms. 
Overall, halo exclusion only results in a small reduction of 
the total ESD on intermediate scales. Due to the arteficial 
sharpness of the break in the 2-halo central term, halo ex- 
clusion introduces a sharp feature in the total ESD at the 
radius corresponding to this break. Although the sharpness 
of this feature is an artefact of our oversimplified treatment 
of halo exclusion, it does not infiuence our overall results. In 
fact, including or excluding halo exclusion has only a small 
impact on the total x^-values of our models. For example, for 
the WMAP3 cosmology, the reduced x^ of our fiducial model 
is 3.1, compared to 4.6 if halo exclusion is ignored. This dif- 
ference is much smaller than that between the WMAPl and 
WMAP3 models. We therefore conclude that our approxi- 
mate treatment of halo exclusion is sufficiently accurate, and 
does not impact our conclusion that the WMAP3 cosmology 
is strongly favored over the WMAPl cosmology. 



r < ri8o(Af). Here M is the average halo mass of the central 
galaxies 
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where ■Pc(M|Li, L2) is the probability that a central galaxy 
with luminosity Li ^ L ^ L2 resides in a halo of mass M, 
and is given by Eq. (Ill|l . The corresponding radius, rigo(M), 
follows from Eq. P^ . 

Although this treatment of halo exclusion is clearly 
oversimplified, we emphasize that previous attempts to in- 
clude halo exclusion in the halo model are also approxima- 
tions (e.g. Magliocchetti & Porciani 2003; Tinker et al. 2005; 
Yoo et al. 2006). In addition, as is evident from Fie:. IA2l halo 
exclusion only has a mild impact on the overall results. The 
black lines, labelled HE, show the ESDs obtained from our 
fiducial model in which halo exclusion is implemented using 
the method outlined above. For comparison, the red lines, 
labelled NOHE, show the results in which we ignore halo ex- 
clusion altogether (i.e. in which the 2-halo terms are simply 
computed using Eqs. [20] and [21]). The dashed lines show 
the corresponding 2-halo central terms, which are clearly 
suppressed on small scales in the HE model. Since brighter 
central galaxies are hosted by more massive (and therefore 
more extended) haloes, the effect of halo exclusion is appar- 
ent out to larger radii for brighter galaxies. Note also that 
the truncation is fairly sharp; this, however, is partially an 
artefact due to our approximate treatment in which we have 
only considered the average halo mass M{L\,L2)- In reality, 
the central galaxies live in haloes that span a range in halo 
masses, and thus a range in sizes. If this were to be taken 
into account, the truncation would still occur at roughly the 
same radius, but be less sharp. 

Although halo exclusion clearly has a strong impact on 
the 2-halo central term, the impact on the total ESD is only 
modest. This mainly owes to the fact that the total signal 
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Figure A2. The ESD is sliown for tliree luminosity bins. The black Unes refer to the fiducial model (HE) and the red lines to the 
model without halo exclusion (NOHE). The solid lines indicate the total signal, whereas the long dashed lines show the 2-halo central 
terms (note that the we ignore halo exclusion for the 2-halo satellite term). Although the 2-halo central term is strongly affected by halo 
exclusion, the impact on the total ESD is only mild. Note that the sharpness of the dip in the black solid linos is (at least partially) an 
artefact of our oversimplified treatment of halo exclusion, as discussed in the text.. 



